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Abstract. In this paper, we present new explicit simultaneous rational approxima- 
tions converging sub-exponentially to the values of Bell polynomials at the points of 
the form (7, l!(2a + 1)C(2), 2!C(3), . . . , (m - l)!(a + 1 + (-lfa)((m)), m = 1,2,..., a, 
a € N. 



1. Introduction 

In 2007, A. I. Aptekarev and his collaborators [HH3] discovered a sequence of rational 
approximations p n /g n converging to Euler's constant 

7 = lim ( 1 + ~ H h - - log fc 

fc->oo y 2 A; 

sub-exponentially. More precisely, the numerators p n and denominators g n of the ap- 
proximations are positive integers generated by the following recurrence relation: 

(16n - 15)g n+ i = (128n 3 + 40n 2 - 82n - 45)g n 

- n 2 (25Qn 3 - 2A0n 2 + QAn - 7)g n _i + n 2 (n - l) 2 (16n + l)g n _ 2 

with the initial conditions 

Po = 0, pi = 2, p 2 = 31, 
g = 1, qi = 3, g 2 = 50 

and having the following asymptotics: 



e V2n 



p~„ - 7 g„ = (2n)!^ f ^ + 0(^; 



( 4e )3/8 

The authors [7] found explicit representations for these sequences: 



(2) ?» = X)u) ^™ = I](j (n + k)\{H n+k + 2H n ^ k - 2H k ) 



fc=0 v 7 fc=0 



ifey words and phrases. Euler constant, zeta value, Bell polynomial, Bernoulli polynomial, Meijer 
G-hmction, rational approximation, saddle-point method. 



here H n = X]fe=i \ ls the n-th harmonic number, H Q := 0. Formulas ([2]) imply that 
q n and p n are integers divisible by n! and 77-, respectively, where D n denotes the least 
common multiple of the numbers 1,2, ... ,n. The linear forms (OQ) do not tend to zero 
even after cancelation of their coefficients by the big common factor and therefore 
this construction does not allow one to prove the irrationality of 7, which is still an 
open problem. Nevertheless, the linear forms ([1]) present good rational approximations 
to Euler's constant 

^-7 = 2vre- 2 ^(l + 0(n- 1 / 2 )) as n -> 00. 

In 2009, T. Rivoal jT/2] found another example of rational approximations to the Euler 
constant 7, viewed as — r'(l), where V is the usual Gamma function. His construction 
is based on the following third-order recurrence: 

(n + 3) 2 (8n + 11) (8n + l% n+3 = (n + 3)(8n + ll)(24n 2 + 145n + 215)y n+2 
- (8n + 27)(24n 3 + 105n 2 + 124n + 25)y n+1 + (n + 2) 2 (8n + 19)(8n + 27)y n , 

which provides two sequences of rational numbers P n and Q n , n > 0, with the initial 
values 

P = -1, Pi = 4, P 2 = 77/4, 

Qo = 1, Qi = 7, g 2 = 65/2 

such that ^ converges to 7. The sequences P n , Q n satisfy the inclusions (see [61 Corol- 
lary 5]) 

n\ Q n , n\ D n P n e Z 
and provide better approximations to 7 

< Coe -9/2nV3 +3 /2nV3 ; |Q n | = O^ 8 "^ 8 ) as n ^ 00. 

Unfortunately, this convergence is not fast enough to imply the irrationality of 7. 

In the same paper [12J, T. Rivoal considered a more general construction which is 
based on simultaneous Pade approximants to Euler's functions 

£ s {z)= —dt, sGN,^C\[0;+oo). 

Jo z — t 

This approach allows one in principle to find rational approximations to some other 
constants related to higher derivatives of the Gamma function r^(l), n > 2. Note 
that T^(l) can be written as (see [31 p. 175]) 

rW(l) = Y n (-j, C(2), -2!C(3), . . . , (-l)»(n - l)\((n)), 

where Y n is the Bell polynomial (see below, for definition and basic properties). Unfor- 
tunately, the details of this construction become rapidly very complicated and T. Rivoal 
presented explicitly rational approximations only for two numbers 7 and ((2) — 7 2 = 
r"(l) — 2r'(l) 2 . More precisely, he constructed a sixth order linear recurrence with 

2 



Pn 
Qn 



polynomial coefficients of degree 25, which has three solutions {ai, n }n>0) { a 2,n}n>o, and 
{b n }n>o such that a\ >n , a 2 ,ni bn ^ 



(3n)!(3n+2)! 



7 



Z and 



C(2)-7 2 - 



«2,r 



< 



16. 



3/8 



/r 



exp(4\/2n 3/4 - 5 



2/8n 



1/4 N 



as n — y oo, 



where Co is some positive constant independent of n. Notice that the better inclusions, 
namely, n\ 2 b n , D n n\ 2 a^ n , D 2 n\ 2 a2, n 6 Z were proved in [6], Corollary 6]. 

Recently, the authors j8] gave a new interpretation of Aptekarev's approximations to 
Euler's constant in terms of Meijer G-functions and hypergeometric-type series. This led 
to finding new rational approximations to 7 generated by a second-order inhomogeneous 
linear recurrence with polynomial coefficients. The denominators and numerators of 
these approximations are given by the formulas 



Qn 



E 



k\(2H n . k -H k 



k=0 v/ k=0 

The sequence {q n } n >o satisfies the second-order homogeneous linear recurrence 



q n+2 - 2(n + 2)q n+1 + (n + l) 2 q r , 







with g = 1) (Zi = 2, and the sequence {p n }n>o is a solution of the second-order inho- 
mogeneous linear recurrence 



with p = 0, pi 



and 



q n+2 ~ 2(™ + 2)q n+1 + (n + l) 2 g n 
1. Moreover, one has 
^- 7 = e- 4 ^(27r + 0(n- 1 / 2 )) 

q n 



n 



q„ 



~TTT [ " 

V: 



+ 0(n 



-1/2^ 



as n — > 00 



as n — y 00. 



nV4 V2v7re 

In this paper, we generalize the construction of work [8] and present explicitly new 
simultaneous rational approximations converging sub-exponentially to the values of 
Bell polynomials at the points of the form (7, l!(2a + l)((2), 2!((3), . . . , (m-l)!(a + l + 
( — l) m a)((m)), m — 1, 2, . . . , a, a £ N. Note that our approach is different from that of 
Rivoal and based on application of Meijer G-functions and complex integrals. 

First we recall several known facts related to the Bell polynomials. The exponential 
(complete) Bell polynomials (first effectively studied by Eric Temple Bell |2] and named 
in his honor) are the polynomials Y n (xi, . . . ,x n ) in an infinite number of variables xi, 
x 2 , ■ ■ ■ , defined by the formal series expansion [3J §3.3]: 



(3) 



cxp 



57 

\m=l 



ml 



E 

ra=0 
3 



Y n (xi, X2, 



, Xr 



ni 



An explicit representation of Y n is given by [21 p. 264] 

Y n (x u ...,*„) = £ • • • Q) kn , 

7r(n) 

where the summation takes place over all partitions 7r(n) of n, i.e., over all n-tuples of 
non-negative integers (ki, . . . , k n ) such that 

n 

J2jkj=n. 

3=1 

The Bell polynomials satisfy the following recurrence relation (see [21 p. 263]): 

Y n+1 (x h . . . ,x n+ x) = ^2 ( l ) x k+iY n -k(xi, ■ • • ,x n - k ), n > 0, Y = l, 

k=0 ^ ' 

which implies immediately that the complete Bell polynomials have integer coefficients 
and therefore, 

n! 

^ ' VMk 1 \2\^k 2 \---n\ k -k n \ E 

for k±, . . . ,k n as above. The first six complete Bell polynomials are as follows: 

Y = 1, Yi(xi) = xi, Y 2 (x 1 ,x 2 ) = x\ + x 2 , Y 3 (x 1 ,x 2 ,x 3 ) = xl + 3xix 2 +x 3 , 

Yi(x 1 , x 2 , x 3 , x 4 ) = x\ + Qxjx 2 + 4xix 3 + 3xj + x 4 , 

Ys(xi,x 2 , x 3 , 24, 25) = x\ + 10x\x 2 + 10x^3 + lhxix\ + 5XiX 4 + 10X223 + X5. 

Let xi, x 2 , . . . , x n , . . . , yi, y 2 , . . . , y n , . . ., be two infinite sequences of independent 
variables. Then from ([3]) we get easily the addition theorem [21 p. 265] 

Y n (xi + Vl , . . . , x n + y n ) = [Y(x) + Y(y)] n . 

In ordinary notation this is equivalent to 

(5) Y n (xi +yi,...,x n + y n ) = ^ ( )Y k (xi, . . . , x k )Y n _ k ( yi , . . .,y n -k)- 

k=0 ^ ' 

The Bell polynomials play an important role in taking the n-th derivative of a composite 
function. Namely, the n-th derivative of the function e^ x ' can be expressed in terms of 
known quantities by 

(6) (4-) e /W = ■ Y n (f(x), f"(x), . . . , /<»>(*)). 



dx / 

This formula is also known as Faa di Bruno's formula for the n-th derivative of the 
composite function. 

Let as usual, (A) m be the Pochhammer symbol (or the shifted factorial) defined by 
(A) = 1, (A) m = A(A + 1) • • • (A + m — 1), m > 1, and H n m ^ be a generalized harmonic 
number given by H n m) = £Li ^ and i# = H n . 



Now we can formulate our main result. 



Theorem 1. Let a > 2 be an integer. For /j, — 1, 2, . . . , a — 1 and any non-negative 
integer n, define the following sequences of rational numbers: 

n / \ a n , \ a 

(7) Qn:=J2{ n k ) kleZ > ^ := E(j k\Y,(n(k),r 2 (k),...,r,(k))eQ, 

h—n \ / h — n \ / 



where 



(8) r m (k) := (m - 1)! (afl£J + (-!)*> - l)^ m) ) , 



k = 0, 1, . . . , n. 



Let 



a M := y M ( Tj l!(2a - 1)C(2), 2!C(3), ...,(/.- l)!(a + (-l)"(a - 1))C(^))- 
Suppose that the coefficients b m (a) are defined by the expansion 
(9) 

" {2-^) m _ N ^ (2- 



-a\og(l+J2 



m=l 



(to + 1)1 



z"' - 



E 

m=l 



a M-l m 

rn! 



^& m (a)z™+0(z a+1 ), 



Id < 1. 



m=l 



In particular, we have &i(a) = —a, &2( a ) = "^p, ^3(0) = ( 1 ~ a ^ a ~ 3 ) . 

Taen /or every /x = 1,2, ... , a — 1 £aere exists a positive constant c M = c^(a) snc/j 
£aat /or any non-negative integer n we have 

'0-1 



|Pn,/i - 5na M | < „ j_ exp > (-1) b m {a) cosl — 
n2 + 2a \ _ — ; v a 



(2ixm\ 1 _ 



n 



m=l 



Moreover, D£ ■ p njfl G Z and t/ie following asymptotic formula holds: 



Qn = — irrr exp ( ^(-1)^(0)^^ ) (1 + 0( n - l /»)) 



as n — )■ 00. 



\m=l 



The sequences {p n ,^/<?n}n>o, /i = 1, 2, . . . , a — 1, provide good simultaneous rational 
approximations converging sub-exponentially to the numbers a M . Note that 

«i = 7, « 2 = 7 2 + (2a - 1)C(2), a 3 = 7 3 + (6a - 3) 7 C(2) + 2C(3), .... 

Corollary 1. Let a > 2 be an integer. Then for /x = l,2,...,a — 1 we have 



< c^exp ^^(-l) m & m (a)(cos(27rm/a) - ljn 1 "^ 



< exp ^a(cos(27r/a) - ljn 1 »(1 + o(l))j , 

where c^ = c^(a) is a positive constant independent of n. 

In particular, for Euler's constant we have 

Corollary 2. Let a > 2 be an integer. Let the sequence {g n } n >o be defined in (GO) and 
Pn = Yl=o (l) a k\(aH n ^ k - (a - l)H k ). Then 



7 



Pr, 



< e 



a(cos(27r/a)-l)ra 1 ~i (l+o(l)) 



5 



Consider several examples. For a = 2, we get rational approximations to Euler's 
constant studied in [8]. 

Example 1. For a = 3, define three sequences 

n / \ 3 n / \ 3 

?» = Eu) ^ = Eu) k\(3H n _ k -2H k ), 

fc=0 ^ ' fc=0 ^ ' 

^ = E u) * ! (^-* - 2H ^ + + 2H ? ] ) ■ 



Then we have 

Pn,l 



7 - 



-9/2n 2 / 3 +3/2n 1 /' i 



7 2 + 5C(2) - ^ 

On 



< c 2 e 



-9/2n 2 / 3 +3/2n 1 '' 3 



and 



n! 



Qn 



„3n 2 /3_ n l/3 / e 



1/3 



5/3 



,2tt\/3 



+ 0(ri 1/3 ) | as n -> oo. 



Applying Zeilberger's algorithm of creative telescoping (see [11] ) it is possible to show 
(by the same way as in [HI Lemma 1]) that the sequences p n>1 , q n are solutions of the 
third-order homogeneous linear recurrence 

(n + l)(8n - 9)/ n+ i = (24n 3 + 13n 2 - 32n - 18) f n 

- n(24n 3 - 75n 2 + 52n - 5)/„_i + n(n - l) 3 (8n - l)/ n _ 2 
with the initial conditions 

Po,i = 0, = 1, p 2 ,i = 13/2, 

<?o = 1, 9i = 2, <? 2 = 11 

and the sequence p„ i2 is a solution of the third-order inhomogeneous linear recurrence 
(n + l)(8n - 9)/ n+ i = (24n 3 + 13n 2 - 32n - 18)/ n - n(24n 3 - 75n 2 + 52n - 5)/ n _i 

+ n(n - l) 3 (8n - l)/ n _ 2 + 2(8n 4 - 17n 3 + 74n 2 - 12n - 9)/(n(n + 1)) 
with the initial values po,2 = 0, pi >2 = 18, p 2)2 = 95. 
Example 2. For a = 4, put 

«» = EU) k] > PM = E(j ^(4^n- fc -3iJ fc ), 

fc=0 ^ ' fe=0 ^ ' 

!(r 2 (A;) +r 2 (fc)), p n , 3 = E ( J + 3ri(*M*) + r 3 (k)), 

where r m (k) is defined in (jSj). Then we have 



P».2 = E 

k=0 



n / \ 4 



7 



Qn 



-4n 3 / 4 +3n 1 / 2 -5/8n 1 /4 



7 2 + 7C(2) - 



Pn,2 



<7r 



< c 2 e 



-4n 3 / 4 +3n 1 /2_ 5 / 8n i/4 



7 3 + 21 7 C(2) + 2C(3)- 



Pn,3 



< c 3 e 



-4n 3 / 4 +3n 1 /2_ 5 / 8n i/4 



and 



C O nl 4n 3 / 4 -3/2r i 1 /2 +5 / 8n i/4 
n 17/8 e 



(l + CKn- 1 / 4 )) 



as n — > 00. 



Applying Zeilberger's algorithm of creative telescoping it is easy to show that the se- 
quences q n , p nt i, p n ,2 satisfy the fourth-order homogeneous linear recurrence 

(n + 2) 2 (729n 4 - 162n 3 - 171n 2 - An + 6)f n+2 

= (2916n 7 + 14661n 6 + 20862n 5 + 947n 4 - 13008n 3 - 2370n 2 + 1320n + 312)/ n+1 

- (4374n 8 - 18468n 7 - 82674n 6 - 85776n 5 - 1 3062n 4 + 24204n 3 + 13528n 2 + 2680n + 1 68) /„ 
+ n 2 (2916n 7 + 28512n 6 + 61848n 5 + 37667n 4 - 12898n 3 - 17463n 2 - 2692n + 398)/ n _i 

- n 2 (n - l) 4 (729n 4 + 2754n 3 + 371 7n 2 + 2084n + 398) / n _ 2 

with the initial conditions 



% 


= 1, 


Qi 


= 2, 


Q2 


= 19, 


?3 


= 250, 


Po,i 


= 0, 


Pi,i 


= 1, 


P2,l 


= 13, 


P3,l 


= 409/3, 


Po,2 


= 0, 


Pl,2 


= 32, 


P2,2 


= 217, 


P3,2 


= 26444/9, 



and the sequence p„ i3 satisfies the fourth-order inhomogeneous linear recurrence 
(n + 2) 2 (729n 4 - 162n 3 - 171n 2 - 4n + 6)/ n+2 

= (2916n 7 + 14661n 6 + 20862n 5 + 947n 4 - 13008n 3 - 2370n 2 + 1320n + 312)/ n+1 

- (4374n 8 - 1 8468n 7 - 82674n 6 - 85776n 5 - 1 3062n 4 + 24204n 3 + 13528n 2 + 2680n + 1 68) /„ 
+ n 2 (2916n 7 + 28512n 6 + 61848n 5 + 37667n 4 - 12898n 3 - 17463n 2 - 2692n + 398)/ n _i 

- n 2 (n - l) 4 (729n 4 + 2754n 3 + 371 7n 2 + 2084n + 398) / n _ 2 

- 6(729n 10 + 2754n 9 - 17424n 8 - 179680n 7 - 490669n 6 - 549106n 5 - 194460n 4 
+ 100424n 3 + 105332n 2 + 30840n + 3184)/(n(n + l) 2 (n + 2)) 

with the initial values p ,3 = 0, pi t 3 = 60, j9 2 ,3 = 402, p 3i3 = 50761/9. 



2. Analytical construction 



Let a > 2 be an integer. Let us consider the function 




n = 0,l,2,..., 



F 




t=k 



7 



Lemma 1. Let a > 2 be an integer. Then for each /x=l,2,...,a — 1 there exist \x 
constants \„ v , v — 1, . . . , fj,, independent of n such that for n = 0, 1, 2, . . . we have 



P 



n,)i 



F n fl = q n and D£ ■ p nfM G Z, where the sequences q n and p n ^ are defined in ([7]). 
Proof. First let us define the function 

f(t) := alogn! - alogT(n+ 1 -t) - (a- l)logT(t+ 1), < t < n. 
Then we have 

(10) f'(t) = a^(n + 1 - t) - (a - l)ip{t + 1) 
and for m > 2, 

(11) / (m) (t) = (-l) mr - 1 aV> (m ~ 1) (ra + 1 - t) - (a - l)^™" 1 ^ + 1). 
Using the well-known formula for the derivatives of i/j(t), 

(12) ij (m \t) = (-l) m+1 m!C(m + l,t), m G N, 
where 

OO 

c(M) = g(JTip 

is the Hurwitz zeta function, we get 

/("0(t) = (-l) m (m - 1)! ((-l) m -V(™, n + 1 - t) - (a - l)((m, t + 1)), m>2. 

We will be interested in the values of f^ m \t) at the integer points t = k, < k < n. 
Then for m = 1 taking into account the well-known properties of the digamma function 

V>(l) = -7, #i + l) = #n-7, «eN, 

by ( TTU|) . we have 

(13) /'(A;) = -7 + aH n „ k - (a - l)H k , 
and from ffTTT) . (fT2j) for any integer m > 2 we get 

(14) /C-)(fc) = (m - l)!(((-l)- 1 (a - 1) - a)C(m) + off^ + (-l) m (a - l)^ m) ). 

Now if denote r m (k) by formula (jSJ) we get D™ ■ r m (k) G Z, and moreover, from f[T3"l) . 
(JIl]) we obtain 



r m(k) 



f'(k) + 7; if m = l; 

/( m )(A;) + (m- l)!(o+ (-l) m (a- l))((m), if m > 2. 



Now let us notice that F(n,t) = e*® and for calculating the /i-th derivative of F(n,t) 
we can apply the Faa di Bruno formula (jH]) to get: 

jX F (n,t) = e /0 • y„(/'(t), /"(*), • • • , /<">(*)), 1 < M < a - 1. 



Now considering the value Y^r^k), r 2 (k), . . . , r^{k)) and applying the addition formula 
we obtain 

Y^n(k),r 2 (k), r M (*)) = £ K (/'(&)> /"(A;), . . . , / W (A;)) 

x V^(7, l!(2a - 1)C(2), 2!C(3), . . . , - 1/ - l)!(a + (-l)^(a - l))C(/i - 1/)). 
This implies that 

r 2 (fc), • • • , r fi{k)) - F M (7, l!(2a - 1)C(2), l)!(a + (-If (a - l))C(/x)) 



x V,( 7 , l!(2a - 1)C(2), ...,(//- 1/ - l)!(a + (-l)^(a - l))C(/x - u)), < k < n. 

Now multiplying both sides of the last equality by e^ k ^ = &!(™) a and summing over 
k — 0, 1, 2, . . . , n we get 

Pn !( u - <7n • ^(7, l!(2a - 1)C(2), . . . , (jjl - l)!(o + (-l)"(a - 1))C(m)) 

= E M V*(7, l!(2a - 1)C(2), -..,(//- 1/ - l)!(o + (-l)^(a - 1))C0* ~ ")) 
i/=i 



x 

or 



where q n and p n ^ are defined in ([7j) and the constants X^ u given by 
K,u = Q ■ W7, l!(2a - 1)C(2), . . . , (// - v - l)!(a + (-l)^(a - 1))C(a* - *0) 

are independent of n. To prove the inclusion Z}^ • p U)/U £ Z, we consider an arbitrary 
monomial of D£ ■ Y" A1 (r 1 (/c), . . . , r m (k)), which has the form 



71 



fci!---ife M ! V 1! / V M ! / l\ kx k x \ ■ ■ ■ /j! fc Mjfe M ! 
x (D n • n(A;)) fcl ■ (D 2 n ■ r 2 (k)) k > ■ ■ ■ (£>£ • r M (A;)) fe ', 

where ki, . . . ,k^ are non-negative integers such that 1 ■ fci + 2 • A; 2 + 3 ■ k 3 + ■ ■ ■ +fi ■ k^ = \i. 
Now by (jlj) and the fact that D™ ■ r m (k) G Z, we get the required inclusion. □ 

Now for each integer < fi < a — 1, and mgK define a complex integral 

/« Ju) := — / F n, t) ( ) e mtu dt 

2m J L \sm7rr/ 



l n./_i' 

(15) 



n + 1, . . . ,n + 1 




where L is a loop beginning and ending at — oo and encircling the points n, n—1, n—2, . . . 
once in the positive direction. Without loss of generality we can assume that L is located 
symmetrically with respect to the real axis. The integral converges according to the 
definition of the Meijer G-function (see [HI §5.2]). Moreover, if /i = a — 1 and \u\ < 1/2 
we can also choose the contour of integration as a vertical line going from c — ioo to 
c + ioo, where c > n is an arbitrary constant. 
Let us also put 



(16) 



n y 
k=0 ^ 



7T \A*+1 



sin7rf 



Antu 



< \i < a- 1. 



Lemma 2. Let a>2, 0</i<a — 1 foe integers. Then 



/n lAt («) + 0(n- a ), 



< /x < a - 2; 
jj = a — 1; 



where the constant in O is absolute 
Proof. First note that 

7T 



F(n,t) 



sin nt 



Airtu 



l) na n\ a . 

a— /t— 1 v ' 



r(t - n) 
r(*+i) 



a-l 



(sin 7rt) 



a— 1 ^iirtu 



Then, since on the segment Ret = — N — 1/2, |Imt| < yo with a sufficiently large integer 
N, we have 



\T(t - n)| < |r(Re (t - n))| = iV - 1/2 -n) 
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T(N + 1/2 + n) 



( e - NlosN+N ), 



r(t-n) 



r(t + i) 



| (sin 7it) 



a— ix— 1 Aivtu I 



0(1), 



it follows that the integral J njAt can be evaluated as a sum of residues at the singular 
points lying inside the loop L. The integrand in (ITS]) has poles of order fi + 1 at the 
points 0, 1, ... ,n and moreover, if // = a — 1, it has additional simple poles at the 
negative integers. Therefore, for fi — a — 1, we have 



J„,«_i(u) = 7„,«_i(u) + (-l) an n! a ^ res 



fc<0 



r(t-n)e 



ITTtu 



t=k \(t(t -!)•••(* - n))"- 1 



/„, a _i(w) + (-l) an n! a ]T res 



fc<0 



~intu 



7T 



j n ,^i(u)+n! a 5i; 



t=fc V r (^ - £ + !)(*(* - 1) • • • (* - n))"" 1 sinvrt 

_-|\(u+l)ft+a-l 



^ (n + £;)!(£;(£; + 1) •••(£; + n))^ 1 



J n , tt _l(u) + 



n + l) a 



k=0 



(_1 U«+l)fe+o-lj^Jo-l 

(^+2)1 

10 



Finally, since 



E 

fc=0 



(n + 2)* 



oo 

^E 



a-l 



x 

fc - IT 



fc=0 v /a fc=0 



e, 



we get the desired assertion. 



□ 



3. Bernoulli polynomials 



The generalized Bernoulli polynomials Bn(x) of order m, where m is a positive 
integer, are defined by the generating formula (see [HI §2.8]) 



(17) 



1)' 



E^'w^ 



\z\ < 2tt. 



n=0 



Numerous properties of these polynomials can be deduced directly from formula (TT7T) . 
We mention here only two of them which will be useful in the sequel. For a detailed 
study of Bernoulli polynomials, see [TU] . For example, comparing powers of z on both 
sides of the equality 



m yz 



z e 



yrn e (y+x)z 



leads to the addition formula 
(18) 



k=0 ^ ' 



Differentiating both sides of (117j) with respect to z and comparing powers of z we get 
the recursion formula 

(19) mB^ +l \x) = (m - n)B^\x) + n{x - m)B ( ™\{x). 
Setting n = m in (fl9|) we get B^n +1 \x) = (x — m)B^}_ 1 (x), which implies 

(20) B^ +1 \x) = {x- l){x - 2) • • • (x - m). 

Lemma 3. 1) Let m be a positive integer. Then the following series expansion: 



sin m z 



E 

n=0 



-\Y^B { Z\m/2) 
{2n)\ 



„2n 



\Z\ < 7T, 



holds. 

2) Let m be a positive even integer. Then 



(21) 



zfTW +i) (—\: =o. 



fc=0 



Proof. Replacing z by 2z in (1171) and using the formula sinh^; = (e z — e z )/2 we get 

Izl < 71. 



sinh m z 



~m„(2x—m)z x ^in 'r, 

— — = Ee»wlf 

n=0 
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Setting x = m/2 in the above equality we have 



(m) (™\ 2 "^' 

sinh m 2 V 2 y n! 

n=0 



Z < 7T. 



Taking into account that zj sinh^; is an even function, we get that for any positive 
integer m, 



(22) 

and therefore, 
(23) 



- D 2n+l ^ 2 



0. 



n = 0,l,2,..., 



(m) /m\ 4 n z 2n 



sinh m z A^-^™ \2J (2n)\ 

n=0 v 7 



\Z\ < 7T. 



Now replacing 2 by in (|23p and recalling that sin 2 = — i sinh(i^), we get the required 
expansion. 

To prove equality (l2~Tj) . we consider the addition formula ( 1181) with m and n replaced 
by m + 1 and m, respectively. Then setting x — 1/2, y — (m + l)/2 we have 



fc=0 



/// \ „( m+ i) / m + 1 



2^ 2 m ^( m +i) 



m + 2 



Now taking into account that m is even and applying ( 1201) . we get5^ +1) (m/2 + l) = 0, 
and the lemma is proved. □ 



4. Properties of the integrals I n ^{u) 



Lemma 4. Let a>2,0</x<a — 1 be integers. Then we have 



4,0(1), 



if fi = 0; 



[M/2] 

J] C 2 j + l^/ n ,2j+l(0), i/ /X 25 Odd; 

5=0 

C2j>Re/ n) 2j(l), i/ > 2 is even; 

where 7^ and c J)At , < j < /x, are real constants independent of n. 
Proof. The proof is by induction on /x. For /i = 0, we deduce easily from (|T6|) : 



n 

7 n , (l) = ^res (F(n,t)e^ 



7T 



fe=0 



sin 7rt 



£*>>*) = £□ fc! = i v 

fc=0 fc=0 ^ 7 



Similarly, for /x = 1, we have 



n ( 2 \ n d 

7„, l( o) = Eg (*"(».«> ' J = S S F( ' M > 

fc=0 v 7 A.-0 



F n A- 



t=k 
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For /x = 2, by ffTBl) . we get 

md = Ejs (*■(».*)«"■ ( s * 

fe=0 v 



sin 7rt 



Then using the expansions (see Lemma [3]) 

(24) f * r' ( n^ ^-'I'M^"^) ,,,,) 

(25) e int = e ink e in(t-k) = ^-(t - k) j + O ((t - k)^ 1 ) , 

3=0 J ' 

(26) F(n, = E 4 F(n, t) ^ + O ((t - k)^) 
in a neighborhood of the integer point t = k we obtain 

and therefore, F n ^ = 2Re/ n 2(l). Now assume that > 2 and the formula holds for 
0, 1, 2, . . . , /i — 1. We will prove it for /i. If /i is odd we have 

n f +i 

fc=0 v 

and from f[2"ll) . f[2"6"j) we deduce that 

1 [/V2] 

Hence, by the induction hypothesis, we conclude that 

[M/2] [M/2]-l 

= /i!/ n , M (0) - At! E d 2j:fl F n:fl _ 2 j = + E c 2m,^n,2m(0), 

i=i z=o 

as required. If /i is even, then by ( !T6|) . we have 

Ul) = Eg %<)^'(^)' ) 

fc=0 v 7 

and using the expansions (!24l) - (l26j) we obtain 

^ (-Ij^vQ^f P(g±l) ^ 
Kev^ij ^ ( 2j -)i 2^ (2/)!(/i-2j-20! n,Ai " 2j " 2 ' 

(27) J "° l ~° 

j'-^n,^ "I - ^ ^ d,2v,iiF n ^—2ui 
^" u=l 
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where di v ^ are some real constants independent of n and 

Taking into account (122]) we can rewrite the last equality as 

which implies g^ iM = 0, by Lemma [H Now according to ff2~Tj) we have 

A*/2— 1 ^ 



1 " /2-1 ~ 

ReJ„ iM (l) = — j-F n>At + d 2v ^F n ^ 



U! 



from which, by the induction hypothesis, the lemma follows. □ 

Lemma 5. Let a, fi be integers satisfying a>2,0</i<a— 1 and u G E. T/ien /or 
eac/i n = 0, 1, 2, . . . we /iave 



4» = 7 ^-^r T £ (" 1 ) fc ( a £ 1 )4, a -i(a + M -/,-l-2A;). 

i — n \ / 



a— /1— 1 

1 

(2tt«" 

Proof. From (1151) and the reflection formula for the gamma function 

(-l) n 7T 

r(* - n)r(n - * + 1) = - — — 

sin Tit 

we have 

~ 2ri J L T^it + l)r a (n-t + 1) VshiTrJ ' / 
:-l) m «! a /" T a (t-n) fsmTTty-"- 1 ^ tUdt 



2m J L r a_1 (t + 1) v 7T 

Now replacing the function sin(7rt) by [e mt — e~ tnt ) / (2i) and opening parenthesis, by 
the binomial theorem, we get 



- {2 ^y, „ 2-j ^ ^ /, ) ], \- HI I 1) 



fc=0 

which completes the proof of the lemma. □ 



Lemma 6. Let a,// 6e integers satisfying a > 2, < fi < a — 1, ana? u G R. T/ien /or 



eac/i n = 0, 1, 2, . . . we have 



where the bar stands for complex conjugation. 
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Proof. Making the change of variable 1 1- > t in the integral I n ^{u) and using the equal- 
ities T(z) = T(z) and sin (2) = sin (2) we get 

1 C f \ - 1 C f \ 

_ ' ' " ITTUt 7T / TT 1 / T\ / " \ i-KUt 



In,M = W~ / F(n, t) — ^ e l ™< dt = -— F(n, t) — e™ ut dt 
2m J L \sm(irt) ) 2m J L \sm(nt) ' 



[ ( y r Y F(n,t)e- W rft = — / F(n,t) f — ^— V +1 e- imi * dt 
2?rz J L \sin(7rt) / V ; 2m J L y 1 VsinTrt/ 



as required. □ 

5. ASYMPTOTICS OF THE INTEGRAL J n>0 _i(u) 

Lemma 7. Let a, it e Z, a > 2, |it| < a, and n be a sufficiently large positive integer. 
Then all a roots of the polynomial 

„ / \ "T«„/ 1 \a a— 1 

p u (r) = e n(r - 1) - r 
are given by the following asymptotic expansions: 

qq m{2itk — iru)i 

/ \ ! , 2 - m a m _i e 
T fc (u = l+> : — — , fc = 0,l,...,a- 1. 

' m! n m/a 

m=l 

Proof. First note that the polynomial p u (r) has not real roots on (— 00, 0]. Indeed, if 
we suppose that r = —x, x > 0, is such a root, then x) = and 

(28) (-l) u n(x + l) a + x a_1 = 0. 

On the other hand, we have that the left-hand side of (12 8j) is positive (negative) if u is 
even (odd), which is a contradiction. Therefore, we can consider the equation 

(29) n{j - l) a e inu = r"" 1 

in the complex r-plane with cut along the ray (—00; 0]. Then it is easily seen that the 
equation ( 129]) is equivalent to a relations of the form 

. . 1 7tu — 2-nk ■ ... . 1 1 

(30) ni e — '(r-^r 1 "", k = 0, 1, . . . , a - 1, 

where r 1- ^ = e( 1 ~^ logT and logarithm is defined by its principal branch. From ( 130|) we 
get immediately all a roots of the polynomial p u (r) : 

(2 Z[ fc = ™- ) _ x exp( 2(27rfc - TO) i) , 
31 r k (u =1+ Pl ° ' + ^+O n- 3 / a , fc = 0,l,...,o-l. 

It is possible to find a complete asymptotic expansion in (l3"Tj) if we apply the Lagrange 
inversion formula (see [U §2.2]) to the equation (|30|) . Indeed, substituting r — 1 = z 
and rewriting (|30|) as 

/qO\ Z - 1 27Tfc-7rU,- 

1 J (z + l) 1 "* 
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we get that there exist positive numbers p\ and P2 such that for n x l a < p\ the equation 
(1321 has just one solution z in the domain \z\ < P2, that is 



(33) z = 



oo 

C m m(27rfc — ttu) 

e a ' 



m=l 

where the coefficients c m are given by the formula 

(34) c m = -((±T~\z + lT^) 

ml \\dz/ j z=Q 

If we suppose that there is another solution of (1321) with \z\ > P2, then since the function 
(|z|+i)i-i/° i ncreases from to +oo as \z\ increases from to +oo, we get 

p* < !f! < !f! = n -v« 

(p 2 + i) 1_1/a ~ (M + i) 1_1/a ~ k + i| 1_1/o 

which is impossible if n is sufficiently large. Hence formulas ( 133]) . ( jMJ) give the complete 
asymptotic expansion for Tk(u), 

T k (u) = i+y. - — — e " ' = 0,l,...,a- 1, 

m=l 

and the lemma is proved. □ 

Lemma 8. Let a, tt G Z, a > 2, |u| < a. T/ien i/ie following formula holds for the 
integral I n>a -i(u) as n — > oo : 

(35) I n , a ^(u) = ( ~ 1)M( ^ ) " g " / ff(T)e»/W dr • (1 + ©(n 1 /"- 1 )), 

in 2 

(36) /(t) = a(r — 1) log(r — 1) — (a — l)r logr — r + r logn + inru, 

(T- 1)5 

(37) <7(t) = 1-^-, 

r 2 

and L\ is a loop beginning and ending at — oo, encircling the points 1, 0, —1, . . . once in 
the positive direction and intersecting the real axis at the point r := To(0) i/u ^ 0, and 
Li is a vertical line going from r — ioo to r + ioo ifu — 0. 

Proof. To make the many- valued functions f(r) and g(r) definite, we consider them in 
the r-plane with cut along the ray (— oo, 1] fixing the branches of logarithms that take 
real values on the interval (1, +oo) of the real axis. 

Note that the contour of integration L in the integral I n ,a-i{ u ) defined in (|15|) is an 
arbitrary loop beginning and ending at — oo and encircling the points n, n — 1, n — 2, . . . 
once in the positive direction or an arbitrary vertical line Re r = c with c > n. Now let 
us denote by r := t (0) = 1 + + 0(n~ 2 / a ) one of the real roots of the polynomial 
Po(t) (see Lemma [7]) and suppose that the contour L (loop or vertical line) intersects 
the real axis at the point utq + 1 and for any point t of L, \t\ > utq + 1. 
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Then the asymptotic expansion of the Gamma function for large \z\ (see [3 §2.11]), 

(38) lo g r(s) = (z - log* - z + X - log(27r) + 0(\z\-\ 

where | argz| < 7r — e, e > and the constant in O is independent of z, implies that 
the following formula holds for the integrand of I n ,a-i{u) on the contour L : 

T a (t-n) M 



e tntu = exp { a i og r(t - n ) - (a - 1) \ogT(t + 1) + mtu} 



r a -!(t + r 

= exp ^a(t — n — log(£ — n) — (a — 1) (t + log(t + 1) + mtu 

+ an-t + a-l +logv / 2^ + 0(7i 1/a_1 )| . 
The change of variable t = nr + 1 yields 



Jirtu _ ( ^u^/^ c an-a,nlagn-(a-3/2)loKn \ T -Q 3 c n f( T )(\ _|_ 0(n l/a ~ 1 )) 



T a (t - n) 



T 2 



o+l 
2 



= ( " 1) ^ ) ' g (r)^)(l + Q(nn), 
n 2 n\ a 

where the functions f(r),g(r) are defined in (l36l) . ( 1371) . This completes the proof of the 
lemma. □ 

The next lemma is devoted to the calculation of the asymptotics of I n ^i{u) using 
the saddle-point method. 

Lemma 9. Let a, u G Z, a > 2, \u\ < a. Then the asymptotic behavior of the integral 
In,a-i{u) as n oo is given by the formula: 

a — 1 Tru(a — l)i / a \ 

/ \ (-l) n "(27r)~ e ^ n \ / r , N -™™ i_™ \ ^, -\i a w 

J„,a-i(«) = - V- — ^exp > 6 m ae ■ n 1 - (1 + 0(n 1/a , 

•where the coefficients b m (a) are rational numbers that can be found explicitly from the 
expansion (J9J). 

Proof. It is sufficient to prove the lemma only for u < 0. The case when u > can be 
reduced to this case by Lemma [6j First we determine the saddle points of the integrand 
( IBT)]) . which are the zeros of the derivative of the function f{j). It is easy to see that 
the zeros of the derivative 

f'(r) = alog(r — 1) — (a — 1) logr + logn + mu 

are simultaneously roots of the polynomial p u (r) defined in Lemma El For the root 
Tk{u), k = 0, 1, . . . , a — 1, by Lemma El we have 

27T fc — 7tu ■ / 2irk — iru ■ 

n(u) - 1 = ^ 1 + izl . LZZ. + ( B -v.) 

Since 

7TU 27rA; — 7TU 7T1i 27T 

< < + 2tt , 

a a a a 
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it is easily seen that r (u) is exactly one zero of the derivative f'(r) among all the roots 
(13TT) . Therefore, tq{u) is the only saddle point of the function e n ^ T \ 

By a similar argument as in the proof of Lemma [6], we can find complete asymptotic 
expansions for the points r satisfying the equation Re/' (t) = 0, which is equivalent to 



(39) , ]Z = n~ 1/a , where z = r - 1. 





z 




\z + l\ 


1-1/a 



Application of the Lagrange inversion formula to the equation 

z 



w 



( z+ l)l-l/a 

shows that there are positive constants <5i and 5 2 such that for \w\ < 5 2 there is only 
one solution z satisfying \z\ < 8\ and this solution is an analytic function of w : 

E(2-m/a) m _i m 
w 

(actually, by the asymptotic formula fl38l) for the Gamma function, the series converges 
if \w\ < a/(a — l) 1_1 / a .) So it is clear that if n > 5%, there is one and only one solution 
of (13"9~j) in the circle |r — 1| < Si given by 

40 r = 1 + V ^ — r> -7T < v? < 7T. 

m=l 

On the other hand, if we suppose that there is another solution with \z\ > 5%, then by 
the same argument as in the proof of Lemma El we get a contradiction if n is sufficiently 
large. So if n is large enough, the continuously differentiable curve 

I f — 1 \ a fi 
Re/'(r) = k>g |T i =0 

represents approximately a small circle with center at r = 1 and radius n~ l l a . The 
points of this curve are given by their complete asymptotic expansion (HOI . Moreover, 
we have Re/'(r) < inside of this curve and Re/'(r) > outside of it. It is easily seen 
that Re/'(Tfc(w)) = for any integers u and k. 

To apply the saddle-point method to evaluation of the integral I n ,a-i(u), we need to 
choose a new contour of integration L 2 passing through the saddle point tq(u). If u = 
we use the original contour Rer = r . Then for r = t + iy, —oo < y < +oo by the 
Cauchy-Riemann conditions, we have 

d d 

— Re/(r + iy) = -Im— /(r + iy) = -aarg(r - 1) + (a - 1) argr. 
dy dr 

Since Re To > 1, for y < we get 



7T 

— — < arg(r — 1) < argr < 



and therefore, 



(41) Re/(r + iy) = (a - l)(argr - arg(r - 1)) - arg(r - 1) > 0. 
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This implies that Re/(r + iy) strictly increases as y increases from — oo to 0. If y > 
then we have 

< arg t < arg(r — 1) < — 

and hence 

j- Re/(r + iy) = (a - l)(argr - arg(r - 1)) - arg(r - 1) < 0. 

Therefore, the function Re/(ro + iy) strictly decreases as y increases from to +oo. 
This proves that Re/(ro + iy) attains its maximum on L\ at the unique point tq and we 
can apply the saddle-point method to calculate the asymptotics of /„ ia -i(0) (see (J48]) 
below) . 

If u 7^ we define L<i as a contour consisting of three parts: 

(i) the half-line r = r + iy, — oo < y < 0; 

(ii) the segment r = r + re lip , if = arg(r (w) — r ), < r < \r (u) — r |, connecting 
the points r and t (u); 

(Hi) the half-line r = — x + ilmTo('u), — Rer (M) < x < +oo. 

If u = —a both parts (ii) and (m) give a ray going from tq to — oo along the upper 
bank of the cut [— oo, 1]. 

Now show that we can replace the contour L\ in the integral I n ,a-i(u) by the contour 
L 2 , i.e., show that 

(42) f e nf{r) g(r)dt = j e n/(r) g(r) dt. 

For this purpose, consider the circle |r| = N, where N is a sufficiently large integer. 
Suppose that L^, (Lf, Lj) are points of intersections of this circle with the contours 
Li, L 2 in the lower half-plane (upper half-plane), respectively. Then to prove (1421 . it 
is sufficient to show that 

(43) / e nf( - T) g(r) dt 0, / e n/(r) #(r) dt -> 

as iV oo, where L^L^ , L\L\ are arcs of the circle of radius iV with center at the 
origin. On the arcs L^L^ and L\L\ of the circle r = Ne %v for iV sufficiently large we 
have the inequalities — 7r < < — 7r/4 and 37r/4 < < ir, respectively (the value (p = n 
corresponds to the upper bank of the cut (— oo, 1]). By Taylor's formula, we obtain 

log(r - 1) = log(iVe^ — 1) = log N + iip — — h 0(N~ 2 ), 

where the constant in 0(N~ 2 ) is absolute. Substituting this expansion in fl36|) . we get 

Re f(Ne iip ) = N log N cos(p)+N log - cos(^) -iV(y?+vra) sm(<p) -a log N-a+0(N- 1 ). 

e 

Note that on the arc L^L^ we have —1 < cos(^) < r / ^/N 2 + Tq , sin(y?) < and 
<y9 + 7tu < 0. This yields 

(44) Re/(r) < r log iV + r log a log iV — a + 0(A^ _1 ) = (r - a) logJV + O(l). 
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Similarly, on the arc L\L\ we have —1 < cos(t^) < — y/2/2, sin(^) > 0, and therefore, 

(45) Re/(r) < -^-NlogN + O(N). 

For the function g(r) on the arcs L^L 2 and L^L 2 , we have the following trivial esti- 
mate: 

(46) \g(r)\ = 0(N 3/2 ~ a ) as iV -)• oo. 

Since the length of each of the arcs L±L 2 and LfL 2 does not exceed 7rN, estimates (14*41) - 
( 1461) imply that the integrals in (pE} are of orders 0(iV ro+1 / 2 - 2a ) and 0(N^ 2 - a ~^^ 2N ), 
respectively Hence, the limiting relations in ( 143]) hold, and we can replace the contour 
of integration L x in the integral I n ^i{u) by the new contour L 2 . 

Now we show that Tq(u) is the unique maximum point of Re/(r) on the contour 
L 2 . Since w < 0, from ( I4T]) it follows that Re/(r + iy) is positive when y < and 
therefore, Re/(r) monotonically increases on the half-line r = tq + £y, — oo < y < 0. 
Similarly, for the half-line r = —a; + % Imro(w) we have 

d d 

— Re/(— x + 2 Imro(w)) = —Re —f(—x + i ImWw)) < 0. 
dx dr 

This shows that the function Re/(r) monotonically decreases on the half-line going 

from To (it) to — oo + i Im To (it). 

Now consider the segment r = r + re tLp defined in {%%). The derivative of Re /(t) on 

this part of L 2 is given by the formula: 

(47) 4- ^ f(r) = Re = Re ( ■ ^] = Re /'(t) cos p - Im /'(t) sin <p. 
dr dr \ dr dr J 

Note that since Re To (it) < Re t , it follows that it/ 2 < <p < n, and therefore the product 
Re/'(T) • cos<£> is positive on the segment. Let us investigate the behavior of Im/' (t) 
on this part of the contour. Note that Im/'(To) = nu < and Im/'(T (u)) = 0. We 
show that Im/'(T) monotonically increases on our segment from ttu to 0. To see this, 
we consider the derivative of Im/'(T). By the Cauchy-Riemann conditions, for r > 0, 
we have 

-f- Im /' (t) = Im (/" (t) • e *) = - Im (/"(r) • (r - t )). 
ar r 

Taking into account that 

f"(r) = *_ a _^l 

T — 1 T 

we obtain 

d 1 fa(r-r ) (a-l)(r-r )\ liar f a(r - 1) (a - 1)t 



Im f(T) = 1 Im f Q(T ~ To) - (a ~ 1)(r - ro) ') 
ar r \ t — 1 r J 



It- II 



Since t is in the upper half-plane, to show that 4- Im /'(t) is positive on the segment, it 
is sufficient to show that the quantity in parenthesis, which we denote by B, is positive. 
For t = T Q + re^, we have 

o(t - 1) (a - 1)t h(r) 



B :-- 

\r — i\ z |t| z |t — 1| z |t| 
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where 



K r ) = ( r o - a)r 2 + 2rr (r - 1) cosy? + r (r - l)(a + r - 1). 

Since r < a, it follows easily that the quadratic polynomial h(r) has two real roots 
r 1; r 2 such that 7*1 < < r 2 and /i(r) is positive on (7*1,7*2) and negative on (— 00,7*1), 
(7*2, +00). Now we show that the point r(u) : = |t (m) — r | belongs to the interval (0, r 2 ). 
Indeed, taking into account that 

\ro ~ 1| = -i- + 0(n- 2 / a ) and |r (u) - 1| = ^ + 0(n' 2 / a ) 
for n sufficiently large, we have \tq(u) — 1| < 2|to — 1| and r < 9/8, and therefore, 

r(u) = \r (u) - tq\ < \t (u) - 1| + |r - 1| < 3(r - 1). 
This implies that 

h(r(u)) > 9(r - a)(r - l) 2 - 6r (r - l) 2 + r (r - l)(a + r - 1) 
= (r - l)(a(9 - 8r ) + 4r (r - 1)) > 

and hence r{u) G (0, r 2 ). This proves that 4- Im /'( r ) i s positive on the segment [0, r(u)] 
and therefore, Im/'(r) increases on it from 7tm to 0, i. e., Im/'(r) < on [0,r(-u)). Now 
by (H71) . we get -pRe/(r) > 011 ^ ne segment connecting r and r (tt) and therefore, 
Re/(r) is monotonically increasing on it. Hence we showed that Tq(u) is the unique 
maximum point on the whole contour of integration L 2 and we can apply the saddle- 
point method (see [5] Ch.6, Th.3.1]) to estimate the integral J nj0 _i(tt). 

Let us recall that when u = —a the contour L 2 consists of the two parts: vertical half- 
line t = r +iy, —00 < y < 0, and horizontal half-line r = — x + iO, x > Rer (a), going 
along the upper bank of the cut (—00, 1]. Let us notice that in this case the function 
g(r) has singularity at the point r = belonging to the contour L 2 . So in order to 
apply the saddle-point method in this case, we should change slightly the contour L 2 to 
avoid this singularity. It is easy to see that instead of the horizontal line r = — x + iO, 
x > —1/2, we can take a contour consisting of the semi-circle r = l/2e l{p , < if < tt, 
and the ray r = — x + iO, x > 1/2. Then we get 

/ e nf{T) g(r)dt= / e n/(r) >g(r) dt + / e nf{T) g(r) dt + 0(C\ 1/2 ), 

J L2 J TQ—ioo J ro+iO 

where C = C(a) is some positive constant independent of n. So we can apply the saddle- 
point method in this case to the contour L 2 a := {r + iy, —00 < y < 0} U {x +i0, 1/2 < 
x < To} and as we can see from fl50|) below the quantity 0{C n n 1 ^ 2 ) does not influence 
on the contribution of the saddle point. 

Finally, applying the saddle-point method we have 
(48) 

= { ~ m ^ e^-t-^"^W)|r(To(7,))]-^(T (7,))e^W)(l + 0(n-i)). 
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In order to find the contribution of the saddle point, we evaluate 

— —)m-l e a 



f(r (u)) = -alog(r (w) - 1) - t q (u) = irui + hgn - 1 - 

m= 

a (0 m+l \ -Krnui 

- alogfl + y 4 + O^- 1 - 1 ^). 

V ^ (m + 1)1 n m / a / v 7 

m=l v y 

Now expanding the logarithm in powers of e~^~ n~ x l a we get 



ml n m / a 

m=l 



V ^ fm + l)! n m / a ) ^ ml 



-< (m + l)l n m / a 
(49) m=1 



m=l 



m=l 

where b m (a) are rational coefficients depending only on a, in particular, bi(a) = —a, 
6 2 (a) = (1 - a)/2, 6 3 (a) = (1 - a) (2a - 3)/(6a). From gSJ it follows that 

(50) e»^("» = IzlC^ exp ( V 6 m (a) e^V^ (1 + 0^)). 

Taking into account that 

g(r (u)) = e^ 1 rT^l + 0^)), f"(r (u)) = ae^ n^(l + O^)), 
we obtain the asymptotic behavior of I n ,a-i{u) '■ 

Vi « = " — V- rrexp h 6 m a e - n 1 » (1 + O n ^ , 



\m=l 

and the lemma is proved. □ 

6. Proof of Theorem 1 

Lemma 10. Let a, \i be positive integers satisfying a > 2, < fi < a — 1. Then there 
exist positive constants Ao = Ao(a), Ai = Ai(a) such that for every positive integer n, 

/ k-i \ , r / m , ^ ^o n! (S, / \ frcm(a-fi- 1)\ x _> 

51 J„ iM < — - -exp > 6 m a cos In 

n2 + 2^ \^ V a / 

\m=l 

and for infinitely many positive integers n the similar lower bound holds: 

'a-l 



Xinl / 7 / \ fumia — u — 1) \ i 



(52) 1/^(0)1 > ^xexp ^6 m (a)cos 

n ~ 2 2a \m=l 

Moreover, the following asymptotic formula 



1/^(1)1 = ^ ^ m ^ exp fE^(«) 008(^^)^-5 ) (l + 0(n-V«)) 
y/a(2Tv) 2 ^n2 + 2a \^ v a / ' 



as n —7- oo takes place. 
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Proof. From Lemma we have 

(53) mi) = ^-^zr a f^-^ (° ~ r j -- i(a - " - 2fc) - 



By Lemmas EH El it follows that the sum on the right of (|53|) contains exactly one term 
with dominant asymptotics, that is I n , a ~i{ a — A*)- Therefore, we have 



p- "' - . + ■ exp f£ 6roW cos(^p^)^-? ] (1 +£>(„-*)), 



as required. Similarly, for the integral In,n(0) by Lemma El we get 

(54) hAO) = ^-^TT if ("!)* (° ~ r J ^- l(a - " - 1 - 2A;) ' 

From LemmaOit follows that if a—fi is odd, then the quantity J„ )/Lt (0) represents a linear 
combination of / n , a -i(0) and sums of complex conjugates I n ,a-i( u ) an d In,a-i{ u ) fo r u = 
2, 4, . . . , a — [L — 1. Similarly, if a — fi is even, then J„ )/Lt (0) is equal to a linear combination 
of differences of complex conjugates I n , a -i(u) and J n)0 _i(w) for w = 1, 3, . . . , a — \i — 1. 
Now by Lemma [9], it is clear that the term with dominant exponent on the right of ( 154|) 
is Re/„ ja _i(a — fi — 1) or Im/ n , i(X _i(a — // — 1) depending on a — /x is odd or even. From 
Lemma [6] we obtain 



(2tt)V n! 

Re J n , a -i(a - fi - 1) = = cos(P(n)) 



x exp ( £ 6 w (a) cos( ^ ~ ^ ~ ) (1 + 0(n^ a )) 

and 



V a 

v.m=l 



a-l 



-l) n (27r)~ n! 
Im J n) a-i(a - /I - 1) = 7= ^— j- sin(P(n)) 

77, 2 + 2a 



x exp ( £ 6 m (a) cosf 7 ^ ~ ^ ~ ) (1 + O^)), 

where 



v,m=l 



7r(a-//-l)(a-l) A. n . /7rm(a - // - 1] 

P(n) = > 6 m (a) sin 

2a ' V a 



-i 211 

n 



m=l 

Since the functions sine and cosine are bounded, we get immediately the required upper 
bound (I5T1) . On the other hand, by Weil's theorem, it is possible to show (see Rivoal's 
argument in the proof of Proposition 13 [12]) that each of the sequences cos(P(n)) and 
sin(P(n)) is dense in the interval [—1, 1]. This implies that there are infinitely many n 
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such that the absolute values of the cosine (sine) are not less than 1/2 and hence the 
lower bound (|52p follows for infinitely many n. □ 



Lemma 11. Let a > 2 be an integer. Then there exists a positive constant A = A (a) 
such that for every \i = 0, 1, . . . , a — 1, we have 

\F n J < ™ exp (^(-lrbM cosily- f) . 

Moreover, the following asymptotic formula holds: 

Qn = — -iL a+ x ex P f E(- 1 ) mfe -( a )^"" J (! + 0{n-^)) 
as n — >■ oo. 

Proof. The proof follows easily from Lemmas HI 121 [TO] and the fact that g n = F nj o = 
4,o(l). ' □ 

Now by Lemmas [TJ [TTJ, the theorem follows. 
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